Search for: All records

Abstract Holsnøy, Norway, offers a world-class natural laboratory for studying the impact of fluid on subducting lower crust. Holsnøy is composed of dry, metastable lower crustal granulite that was infiltrated by fluids along shear zones and seismic fractures during subduction. The infiltration facilitated the localized growth of eclogite facies mineral assemblages along the fluid flow pathways. The duration of the eclogite facies metamorphism, however, remains uncertain. Previous garnet diffusion chronometry studies have estimated timescales ranging from hundreds of years to millions of years based on diffusional relaxation between metastable granulite facies garnet cores and eclogite facies garnet rims and fractures. The shorter timescales are inferred from extremely sharp Ca gradients across chemical contacts present in some garnets whereas the longer timescales are from wider Mg and Fe profiles present in all garnets. The different timescale estimates have led to divergent models for the region’s tectonometamorphic evolution. Here we show that the sharp Ca contacts can be explained by diffusion-induced compositional stress. As Ca is significantly larger than Mg and Fe, its movement strains the crystal lattice and generates stress that limits the relaxation of sharp chemical contacts. When compositional stress is accounted for, the sharp contacts yield timescales that are consistent with the wider Mg and Fe diffusion profiles. We determine that eclogite facies conditions (670–700 °C, 1.5–2.2 GPa) lasted a maximum of c. 300 kyr. The relatively short duration of eclogite facies conditions requires that multiple transient heating events were superimposed on a longer (>106 yr) overall timescale of metamorphism. Granulite facies garnet cores are surrounded by multiple generations of eclogite facies rims formed by interface-coupled dissolution–reprecipitation (ICDR) reactions. The garnet rims indicate two rapid, regional-scale fluid pulses and additional smaller, more localized pulses. The fluid pulses may be linked to episodes of seismic moment release as well as transient heating via exothermic hydration reactions and/or shear deformation. Our model results predict up to 400 MPa of differential stress at the garnet core–rim contacts, consistent with observed eclogite facies microfractures that extend into relic granulite facies garnet cores. The microfractures indicate that ICDR was aided by compositional stress: diffusion ahead of the reaction front generated stress and fracturing that created porosity for further ICDR. Thus, compositional stress can markedly impact both diffusion systematics and intracrystalline deformation. Together, these results show that despite their brevity, transient thermal, fluid flux, and/or baric episodes may exert the primary controls on the mineralogical and rheological development of subducted lithologies, in contrast to the long, slow burial and exhumation typically envisioned for regional metamorphism. 

Abstract Understanding the behavior of chalcophile elements during the evolution of arc magmas is critical to refining models for the formation and distribution of porphyry copper deposits used in mineral exploration. Because magmas in continental arcs undergo copper depletion during their early differentiation, a widely held hypothesis posits that the removed copper is locked at the base of the crust in copper-rich cumulates that form due to early sulfide saturation. Testing this hypothesis requires direct evidence for such copper-rich reservoirs and a comprehensive understanding of the mechanisms driving sulfide saturation. Interaction between oxidized magmas and reducing crustal material in island arcs has been shown to be an efficient process causing sulfide saturation. However, the extent to which crustal assimilation impacts the flux of chalcophile elements during magmatism in thick continental arcs remains to be established. Here, we provide a deep perspective into these problems by studying a suite of subarc cumulate rocks from the Acadian orogen, New England (USA). These cumulates record the imprint of subduction zone magmatism and represent the residues left behind during the genesis of intermediate to evolved Acadian magmas (ca. 410 Ma). We find that the most primitive Acadian cumulates are enriched in copper (up to ~730 µg g–1) hosted by sulfide phases, providing direct evidence for the formation of lower crustal copper-rich reservoirs. The Acadian cumulates reveal a wide range of δ34S values, from –4.9‰ in the ultramafic rocks to 8‰ in the most evolved mafic rocks. The negative δ34S values observed in the most primitive and copper-rich cumulates (avg –3‰) reflect the assimilation of isotopically light sulfur from surrounding sulfidic and graphite-bearing metasedimentary rocks (δ34S of –19 to –12‰), whereas the more evolved cumulates with positive δ34S signatures may have formed from different magma batches that experienced less sediment assimilation. The assimilation of these reducing metasedimentary rocks caused a critical drop in oxygen fugacity (~DFMQ –2.5 to –1.9; FMQ = fayalite-quartz-magnetite buffer) in the evolving magmas, ultimately leading to extensive sulfide saturation and the consequent formation of copper-rich subarc cumulates. Assimilation-driven sulfide saturation may be a common process at the root of thickened arc crusts that triggers the formation of lower crustal copper-rich reservoirs, which play a pivotal role in the fate of copper during arc magmatism. Thus, deeply buried reducing metasedimentary crustal material at the base of continental arcs can act as a barrier to the magmatic flux of chalcophile elements and may play a crucial role in the genesis and distribution of porphyry copper deposits. 

Abstract Gibbs free energy, the fundamental thermodynamic potential used to calculate equilibrium mineral assemblages in geological systems, does not apply to non‐hydrostatically stressed solids. Consequently, there is debate over the significance of non‐hydrostatic stress in petrological and geophysical processes. To help resolve this debate, we consider the effects of non‐hydrostatic stress on the polymorph pairs kyanite/sillimanite, graphite/diamond, calcite/aragonite, and quartz/coesite. While these polymorphs are most relevant to metamorphic processes, the concepts developed are applicable to any single‐component solid reaction. We quantitatively show how stress variations normal to an interface alter equilibrium temperatures of polymorph pairs by approximately two orders of magnitude more than stress variations parallel to an interface. Thus, normal stress controls polymorph stability to first order. High‐pressure polymorphs are expected to preferentially nucleate normal to and grow parallel to the maximum stress and low‐pressure polymorphs, the minimum stress. Nonetheless, stress variations parallel to an interface allow for the surprising possibility that a high‐pressure polymorph can become more stable relative to a low‐pressure polymorph as stress decreases. The effects of non‐hydrostatic stress on mineral equilibrium are unlikely to be observed in systems with interconnected, fluid‐filled porosity, as fluid‐mediated reactions yield mineral assemblages at approximately constant pressures. In dry systems, however, reactions can occur directly between elastic solids, facilitating the direct application of non‐hydrostatic thermodynamics. Non‐hydrostatic stress is likely to be important to the evolution of metamorphic systems, as preferential orientations of polymorphic reactions can generate seismicity and may influence fundamental rock properties such as porosity and seismic anisotropy. 

Abstract The fate of subducted CO₂remains the subject of widespread disagreement, with different models predicting either wholesale (up to 99%) decarbonation of the subducting slab or extremely limited carbon loss and, consequently, massive deep subduction of CO₂. The fluid history of subducted rocks lies at the heart of this debate: rocks that experience significant infiltration by a water-bearing fluid may release orders of magnitude more CO₂than rocks that are metamorphosed in a closed chemical system. Numerical models make a wide range of predictions regarding water mobility, and further progress has been limited by a lack of direct observations. Here we present a comprehensive field-based study of decarbonation efficiency in a subducting slab (Cyclades, Greece), and show that ~40% to ~65% of the CO₂in subducting crust is released via metamorphic decarbonation reactions at forearc depths. This result precludes extensive deep subduction of most CO₂and suggests that the mantle has become more depleted in carbon over geologic time. 

Abstract Inclusions of relic high-pressure melts provide crucial information on the fate of crustal rocks in the deep roots of orogens during collision and crustal thickening, including at extreme temperature conditions exceeding 1000 °C. However, discoveries of high-pressure melt inclusions are still a relative rarity among case studies of inclusions in metamorphic minerals. Here we present the results of experimental and microchemical investigations of nanogranitoids in garnets from the felsic granulites of the Central Maine Terrane (Connecticut, U.S.A.). Their successful experimental re-homogenization at ~2 GPa confirms that they originally were trapped portions of deep melts and makes them the first direct evidence of high pressure during peak metamorphism and melting for these felsic granulites. The trapped melt has a hydrous, granitic, and peraluminous character typical of crustal melts from metapelites. This melt is higher in mafic components (FeO and MgO) than most of the nanogranitoids investigated previously, likely the result of the extreme melting temperatures—well above 1000 °C. This is the first natural evidence of the positive correlation between temperature and mafic character of the melt; a trend previously supported only by experimental evidence. Moreover, it poses a severe caveat against the common assumption that partial melts from metasediments at depth are always leucogranitic in composition. NanoSIMS measurement on re-homogenized inclusions show significant amounts of CO2, Cl, and F. Halogen abundance in the melt is considered to be a proxy for the presence of brines (strongly saline fluids) at depth. Brines are known to shift the melting temperatures of the system toward higher values and may have been responsible for delaying melt production via biotite dehydration melting until these rocks reached extreme temperatures of more than 1000 °C, rather than 800–850 °C as commonly observed for these reactions. 

Abstract Metamorphic devolatilization of subducted slabs generates aqueous fluids that ascend into the mantle wedge, driving the partial melting that produces arc magmas. These magmas have oxygen fugacities some 10–1,000 times higher than magmas generated at mid-ocean ridges. Whether this oxidized magmatic character is imparted by slab fluids or is acquired during ascent and interaction with the surrounding mantle or crust is debated. Here we study the petrology of metasedimentary rocks from two Tertiary Aegean subduction complexes in combination with reactive transport modelling to investigate the oxidative potential of the sedimentary rocks that cover slabs. We find that the metasedimentary rocks preserve evidence for fluid-mediated redox reactions and could be highly oxidized. Furthermore, the modelling demonstrates that layers of these oxidized rocks less than about 200 m thick have the capacity to oxidize the ascending slab dehydration flux via redox reactions that remove H₂, CH₄and/or H₂S from the fluids. These fluids can then oxidize the overlying mantle wedge at rates comparable to arc magma generation rates, primarily via reactions involving sulfur species. Oxidized metasedimentary rocks need not generate large amounts of fluid themselves but could instead oxidize slab dehydration fluids ascending through them. Proposed Phanerozoic increases in arc magma oxygen fugacity may reflect the recycling of oxidative weathering products following Neoproterozoic–Palaeozoic marine and atmospheric oxygenation.